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We construct integrals of motion for multidimensional classical systems from 
ladder operators of one-dimensional systems. This method can be used to 
obtain new systems with higher order integrals. We show how these integrals 
generate a polynomial Poisson algebra. We consider a one-dimensional sys- 
tem with third order ladders operators and found a family of superintegrable 
systems with higher order integrals of motion. We obtain also the polynomial 
algebra generated by these integrals. We calculate numerically the trajecto- 
ries and show that all bounded trajectories are closed. 



1 Introduction 

Over the years, many articles were devoted to superintegrability [1-15]. 
For a review of superintegrability of two-dimensional systems we refer the rea- 
der to the Ref.14. The relation between constants of motion and ladder opera- 
tors in classical and quantum mechanics was acknowledged by several authors 
[3,5,6,16-24]. Detailed discussions of the relation between integrals and ladder 
operators for the two-dimensional harmonic oscillator, anisotropic harmonic 
oscillator and Kepler- Coulomb systems were done [17,18]. Ladder operators 
are more used in context of quantum mechanics. They can provide the wave 
functions and the energy spectrum of the corresponding Schrodinger equa- 
tion and the eigenstates of the annihilation operator are related to coherent 
states [25]. In quantum mechanics, these raising and lowering operators are 
also related to supercharges and supersymmetric quantum mechanics. The 
quantum superintegrable systems with third order integrals of motion [15,16] 
were related to supersymmetric quantum mechanics (SUSYQM) [26,27] and 
higher order supersymmetric quantum mechanics (HSQM) [28,29,30]. A me- 
thod to generate quantum superintegrable systems from supersymmetry was 
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presented in Ref.31. This method allows to generate systems with higher or- 
der integrals of motion. In a recent article, we discussed how ladder operators 
can be used to generate higher order integrals of motion and superintegrable 
systems in context of quantum mechanics [32]. We imposed the separation 
of variables in Cartesian coordinates and the order of the ladder operators 
were arbitrary. These relations between quantum superintegrable systems, 
integrals of motion, polynomial algebras, ladders operators and supersym- 
metry are interesting and provide new insight. In the light of these results 
the study of systems with ladder operators appears to be important also in 
regard of superintegrable systems. The classification of systems with first or 
second order ladder operators in E 2 was discussed [33]. Systems with third 
[16,29] and also fourth [30] order ladder operators were discussed in context 
of supersymmetric quantum mechanics. 

The purpose of this paper is to discuss how the method developped in context 
of quantum mechanics to obtain integrals of motion and polynomial algebras 
from ladder operators [32] can be applied in classical mechanics. The method 
allows to obtain multidimensional superintegrable systems however we will 
focus on two-dimensional systems. We will also discuss systems with third 
order ladder operators. To our knowledge, the study of classical systems with 
higher order ladder operators is also an unexplored subject. We will point out 
that ladder operators appear important in regard of classical superintegrable 
systems. 

Let us present the organization of this paper. In Section 2, we will present a 
method to generate higher order integrals of motion and new classical super- 
integrable systems from one-dimensional systems with ladder operators. We 
present the general polynomial Poisson algebra obtained from these integrals 
of motion. In Section 3, we consider a system that we studied in an earlier 
article concerning classical superintegrable systems in two-dimensional Eu- 
clidean space separable in Cartesian coordinates with a second and a third 
order integrals [17]. We show that this system possesses third order ladder 
operators. We use these operators and results of Section 2 to generate new 
superintegrable systems. We present their integrals of motion and polynomial 
Poisson algebras. We obtain their trajectories. They are deformed Lissajous' 
figure [34]. We show also that all bounded trajectories are closed [35]. 
Before proceeding with the results let us recall a few important definitions. In 
classical mechanics a Hamiltonian system with Hamiltonian H and integrals 
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of motion X a 

H = ^gikPiPk + V(x,p), X a = f a (x,p), a=l,...,n-l , (1.1) 

is called completely integrable (or Liouville integrable) if it allows n inte- 
grals of motion (including the Hamiltonian) that are well defined functions 
on phase space, are in involution {H,X a } p = 0, {X a ,X&} p = 0, a,b=l,...,n- 

1 and are functionally independent ({, } p is a Poisson bracket). A system 
is superintegrable if it is integrable and allows further integrals of motion 
Yb(x,p), {H,Yb} p = 0, b=n,n+l,...,n+k that are also well defined functions 
on phase space and the integrals {H, X\, X n -i, Y n , Y n+ j.} are functio- 
nally independent. A system is maximally superintegrable if the set contains 
2n-l functions. The integrals Yj, are not required to be in evolution with 
Xi,...X n -i, nor with each other. 

2 Ladder operators and integrals of motion 

Let us consider a classical two-dimensional Hamiltonian separable in Car- 
tesian coordinates 

H(x 1 ,x 2 ,P 1 ,P 2 ) = H 1 (x 1 ,P 1 ) + H 2 (x 2 ,P 2 ), (2.1) 

for which polynomial ladder operators (A Xi and A£.) exist. These operators 
satisfy the relations 

{ffi,<}, = ^> {Hi,A-} p = -\ Xi A- , (2.2) 

{A^AtJ^PM), A~.A+. = A+.A~ t = Qi(Hi), < = 1,2 (2.3) 

where Pi(Hi) and Qi(Hi) are polynomials. These relations are the classical 
analog of relation imposed in Ref.32. They are satisfy for many well known 
superintegrable systems such the Harmonic oscillator and the Smorodinsky- 
Winternitz potentials that allow separation of variables in Cartesian coor- 
dinates. From the relation (2.2) the operators fi = A^^A'™ 2 and f 2 = 
A'^A^™ 2 Poisson commute with the Hamiltonian H given by Eq.(2.1) if 
mi\ Xl — rri2^x 2 — with mi,m 2 G Z + . The following sums are also polyno- 
mial integrals of the Hamiltonian H 

h = A%*A£*-AZ*A£*> h = A£*AZ» + AZ*A£* ■ (2.4) 
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By construction the Hamiltonian has the following second order integral from 
separation of variables. 

K = H X - H 2 . (2.5) 

The Hamiltonian H given by Eq.(2.1) is thus superintegrable. 
We will now interested by the algebraic structure generated by these inte- 
grals. In quantum mechanics quadratic [36,37], cubic [14-16] and higher order 
polynomial algebras [32] can be written as deformed oscillator algebras [38]. 
The Fock type unitary representations can be used to obtain the energy spec- 
trum. This is the classical equivalent of the polynomial algebra obtained in 
Ref. 24. The Eq.(2.2) and (2.3) are classical analog of deformed oscillator 
algebra [38]. Such algebras were discussed by A.V.Tsiganov in Ref. 39. The 
Eq(2.3) can also be interpreted as classical analog of the factorization me- 
thod in supersymmetric quantum mechanics [40]. We construct from integrals 
given by Eq.(2.4) and (2.5) the following polynomial Poisson algebra 

{K, h} p = 2A/ 2 , {K, I 2 } p = 2Xh, {h, I 2 } p = 2Q 1 (^(H+K)) m ^ 1 (2.6) 
Q 2 ( 1 -(H-K)r^[mlQ/-(H+K))p/^ 

This is the classical analog of the algebra obtained in Ref. 32. We can relate 
the polynomial by the following relations. The polynomial algebra of su- 
perintegrable systems in classical mechanics play an important role in their 
classification [41]. 

The method can also be extended in N-dimensions by forming the following 
integrals 

Iij = A+p A~p - A~™< A+p , Jij = A + p A~p + A~p A+p , (2.7) 
K ij = H Xi -H Xj , 0<i<j<N. 

3 Construction of new superintegrable systems 

Let us present a system obtained in Ref. 12 and studied in Ref. 13. 

p2 p2 , ,2 

H = ^ + ^ + ^-x 2 2 + V(x 1 ). (3.1) 
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where the potential V(xi) satisfies a quartic equation 

(3.2) 



cj 4 3oj 6 
9^(xi) 4 + 14w 2 x?1/(xi) 3 + (6d- 15— x 4 )\/(xi) 2 + (— xf - 2du 2 x 2 1 )V(x 1 ) 



+ (cxl -d 2 - d^-xt - ^x?) = 0. 
z lb 

This Hamiltonian has two integrals 

p2 p2 , ,2 

^ = ^-^-y^ + ^i), (3.3) 

B = -x 2 P 3 + x 1 P 2 P 2 + (^x 2 1 -3V(x 1 ))x 2 P 1 -^(^-xj-3V(x 1 ))V Xl (x 1 )P 2 . 

z ui z 

In the quantum case V satifies a fourth order differential equation [16] 

nV (4) (x!) = VitfxxV'ixi) + 6(y 2 (x 1 ))" - 2o; 2 x 2 \/ /, (x 1 ) + 2u; 4 x 2 (3.4) 

that can be solved in terms of the fourth Painleve transcendent [42] . In gene- 
ral, Eq.(3.2) has 4 roots and the expressions for them are quite complicated. 
A special case occurs if o> 2 ,c and d satisfy c = 23 "g 63 and d = <i ^- where b is 
an arbitrary constant. Then Eq.(3.2) has a double root and we obtain 

u 2 _ , 



V( Xl ) = — {2b + 5xi + e4x 1 ^/b + xl) (3.5) 

The potentials V(xi) is a deformed harmonic oscillators. The Hamiltonian 
given by Eq.(3.1) with potential given by Eq.(3.5) reduces to isotropic har- 
monic oscillator or anisotropic harmonic oscillator (with ratio 3 :1) when 
6 = 0. For V(x\) satisfying Eq.(3.5) the cubic algebra is 

{A,B} P = C, {A, C} p = -4u 2 B, {B, C} p = 8A 3 + 12HA 2 (3.6) 
-4H 3 - 4 A + 
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3.1 Deformed isotropic and anisotropic harmonic oscil- 
lator 

Let us consider the following one-dimensional systems with potential gi- 
ven by Eq.(3.5) 
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H x = y + ^(2&i + 5xj + et^/h + x 2 ). (3.7) 



. Like its quantum equivalent, this system has third order creation and an- 
nihilation operators 

< = - i ux 1 P 2 + (^f + ^ + Ic^y^T^)^ (3.8) 

+i(-^?u; 3 xi - |L 3 x 3 - ^b\e lS Jb\ + x\- ^e^xly/b 2 + x\) 

a" = P 3 + zu^P 2 + + ^ + ^eia; 2 xi^6f + (3.9) 

-i(--%u 3 Xl - ^c/z 3 - ^bj + xj- ^e^xl y^ 2 + x\) . 

satifying 

{H,a+ i } p = iua+ i , {#,a~} p = -twa". (3.10) 

The results of Section 2 allows us to form the following N-dimensional clas- 
sical superintegrable system 

N P 2 u 2 k 2 i 

It has ladder operators of the same form given by Eq.(3.8) and (3.9) in each 
axis. The polynomial Pj(Hj) and Qj(Hj) are given by 

16 2 
Pj{Hj) = 2Aiuk j H 2 +—i(-b ] u 3 k 3 +b 2 u 3 k 3 )H J + — t(Ab 2 u 5 k 5 J -8b 3 u 5 k 5 J +3b 4 j u 5 k 5 j ) 

(3.12) 

QAHj) = ^(18#; + (6i - 2)b j <Jj) 3 {9H j + 6,(26, - l)o, 2 ) (3.13) 

For the case N=2 the integrals are thus given by Eq.(2.4) and (2.5). These 
integrals generate the polynomial Poisson algebra given by Eq.(2.7) with 
P 1 (H 1 ),P 2 (H 2 ), Q^Hx) and Q 2 (H 2 ) given by Eq.(3.12) and (3.13). The in- 
tegrals Ii and I 2 are polynomials in the momenta of order 3 m i+ m 2 _ \ an d 
3 mi+m2 respectively. The order of the polynomial algebra is 2 x 3 m i+ m 2-i_ 
In the N-dimensional case the integrals are given by Eq.(2.7). Only 2N-1 of 
these integrals are functionally independent and this system is maximally su- 
perintegrable. The integrals 1^-, Jjj and Kij are polynomials in the momenta 
of order 3 m »+ m J — 1, 3 m »+ m j an( j 2 respectively. 
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Figure 1 



Figure 2 



Figure 3 



Figure 4 



The trajectories are obtained numerically directly from the equations of mo- 
tion. We present trajectories for the cases N=2 and N=3 for specific para- 
meters (Figl, Fig2, Fig3 and Fig4). The bounded trajectories are closed and 
correspond to deformed Lissajous's figures. 



4 Conclusion 

Other methods to obtain superintegrable systems were discussed recently : 
the method of coupling constant metamorphosis in context of higher order 
integrals of motion [43] and the method of symmetry reduction [44,45]. In 
this paper, we constructed for a class of classical systems integrals of motion 
from ladder operators. This method allows to generate new classical superin- 
tegrable systems with higher order integrals of motion from one- dimensional 
Hamiltonian for which we know the ladder operators. From the requirement 
that ladder operators satisfy the classical analog of deformed oscillator al- 
gebras we obtained the polynomial Poisson algebra for the superintegrable 
systems. 

We considered the one-dimensional system given by Eq.(3.5) for which the 
potential satisfy a quartic equation. This system has third order ladder ope- 
rators. From this systems and the ladder operators and the method of Section 
2, we constructed the integrals of motion and polynomial Poisson algebra for 
a new two-dimensional suprintegrable systems. A family of superintegrable 
systems in N dimensions can be generated. We present also the trajectories 
in the two and three dimensional case and obtain deformed Lissajous figures. 
All bounded trajectories are closed. 



7 



We discuss in this paper systems with third order ladder operators. Classical 
systems with higher order ladder operators were not systematically studied. 
The study of such systems could provide new superintegrable systems with 
higher integral of motion. 
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Figure captions 

Figl.A trajectory for V = ^(2&i + 5a: 2 + e^xVh + x 2 ) + ^(2b 2 + 5y 2 + 
e 2 Ay^/b 2 + V 2 )- Parameter e x = 1, e 2 = 1, a; = 3, fci = l,fc 2 = 3,6i = 3, 
h = 5, v xo = 1, x = 1, w yo = -3, y a = 1, t=[0,20]. 

Fig2.A trajectory for V = ^f{2bi + 5x 2 + ei 4xV6i + x 2 ) + ^(2& 2 + 5y 2 + 
^V\/b 2 + y 2 )- Parameter e\ = 1, e 2 = 1, u = 3, ki = 3,k 2 = 4,&i = 3, 
&2 = 5, v xo = 1, x Q = 1, = -3, y D = 1, t=[0,20]. 

Fig3.A trajectory for V = ^{2b x + 5x 2 + e^xV^i + x 2 ) + ^(2& 2 + 5y 2 + 

^yVh + y 2 ) + ^(26 3 + 5z 2 + e 3 4zVb 3 + z 2 ). Parameter e x = 1, e 2 = 1, 
e 3 = = 3, h = 7, A; 2 = 11, k 3 = 4, b x = 3, 6 2 = 5, b 3 = 7, w xo = 1, x Q = 1, 
V = -3, y G = 1, Zo = 1, u* = 2, t=[0,20]. 

Fig4.A trajectory for V = ^(2&i + 5x 2 + e^xV^i + x 2 ) + 4?( 2f, 2 + 5?/ 2 + 

e 2 4:y^/b 2 ~+y 2 ) + ^^(263 + 5z 2 + e 3 Az^/b 3 ~~\ r ^ 2 ) ■ Parameter e x = 1, e 2 = 1, 
e 3 = 1, w = 3, ki = 5, fc 2 = 6, fc 3 = 2 61 = 3, b 2 = 5, b 3 = 7, v xo = 1, x = 1, 
%> = -3, y = 1, z = 1, w 20 = 2, t=[0,20]. 
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